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On non-coplanar Hohmann Transfer using angles as 
parameters 

Angel Rincon^’^ • Patricio Rojo^ • Elvis Lacruz^ • 

Gabriel Abelian^ • Sttiwuer Diaz® 


Abstract We study a more complex case of Hohmann 
orbital transfer of a satellite by considering non- 
coplanar and elliptical orbits, instead of planar and 
circular orbits. We use as parameter the angle between 
the initial and transference planes that minimizes the 
energy, and therefore the fuel of a satellite, through 
the application of two non-tangential impulses for all 
possible cases. We found an analytical expression that 
minimizes the energy for each configuration. Some rea¬ 
sonable physical constraints are used: we apply im¬ 
pulses at perigee or apogee of the orbit, we consider 
the duration of the impulse to be short compared to 
the duration of the trip, we take the nodal line of three 
orbits to be coincident and the three semimajor axes 
to lie in the same plane. We study the only four possi¬ 
ble cases but assuming non-coplanar elliptic orbits. In 
addition, we validate our method through a numerical 
solution obtained by using some of the actual orbital el¬ 
ements of Sputnik I and Vanguard I satellites. For these 
orbits, we found that the most fuel-efficient transfer is 
obtained by applying the initial impulse at apocenter 
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and keeping the transfer orbit aligned with the initial 
orbit. 
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1 Introduction 


In 1925, Hohmann studied the transfer between copla- 
nar circular orbits and found that the minimum fuel 
transfer in a Newtonian gravitational field occurs when 
two impulses are applied producing an elliptic transfer 
orbit which is tangent to both of the terminal circu¬ 
lar orbits. A first impulse is used to set the vehicle 
into the elliptic transfer orbit, while a seco nd impulse 
leads t o a circular orbi t at the final radius ( Hohmam] 
( 196flli ; Prussinj ( 1992 )). Many researchers have made 
contributions to the improvement and understanding of 
this type of orbital transfers. More recently, Hohmann 
transfer has been g eneralized from the origina l idea to 
more general cases: Broucke and Prad3 (1994) consid¬ 
ered 7V-impulse transfers between any two coplanar or¬ 
bits (VfV ^4) and for the two-impulse maneuver de¬ 
veloped optimality conditions that lead to a non-linear 
system of three equations and three unknowns, whereas 
Arlulkar and NaiH ( 2012 ) discussed Hohmann transfer 
between two circular orbits but including a dynamical 
approach Lambert solution (i.e., considering the trans¬ 
fer time of the or bit to change f r om on e point to an¬ 
other). Reference Mabsout et al. ( 20091) addressed the 
optimization of the orbital Hohmann transfer consid¬ 
ering only the coplanar case using as optimization pa¬ 
rameter the eccentricity of the transfer orbit. On the 
other hand, different techniques of standard optimiza¬ 
tion had been used for minimizing a cost function. For 
example, one studied the isoperimetric problem of find¬ 
ing the extremal transfers for the given characteristic 
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velocity for the orbits (|KirDichn ikov et al.ll2003i) . Some 
of the authors of this paper ( Lacrud 2010 ) have cal¬ 
culated the solution for the generalized non-coplanar 
Hohmann transfer only for the first configuration (there 
exist four configurations that minimiz e the energy ac¬ 
cording to iKamel and SolimanI (|l999l) l. In this paper 
we consider elliptic orbits and N = 2 impulse transfers 
( Broucke and Pradol 1994 1: taking a split between ini¬ 
tial and transference plane s (b:t 0°), th i s imp roves 
and generali zes the work of Mabsout et al. ( 2009ll and 
differs from iKiroichnikov et al.l (1200311 since we will 
not consider the launch time of spacecraft. We con- 
si der non-coplanar orbits, w hich improve the solution 
of lArlulkar and NaikI (|2012ll . It is relevant to discuss 
the role of orbital transfer in astronomy and engineer¬ 
ing: the orbital transfers are required for a standard 
space mission. Generalized coplanar Hohmann transfer 
had been used to model a space vehicle traveling in el- 
liptic orbits of the Earth and Jupiter around the Sun 
(|Kamel et al.l 1201 111 and shows the importance of this 
kind of study in astronomy and the planetary sciences. 
The standard (non-perturbed) transfer between orbits 
is treated using Kepler problem theory, i.e. consider¬ 
ing Keplerian orbits, because these are non-perturbed 
solutions of the two body problem as a first approx¬ 
imation to a typical orbital motion. We choose to 
work with the standard set of inertial orbital elements: 
O := {a, e, i, H, w, r}, where a is the semimajor axis, e 
is the eccentricity, i is the inclination, H is the longi¬ 
tude of the ascending node and w is the argument of 
periapsis. Finally, the sixth parameter is the epoch t 
indicating the time at which the orbiter passes through 
periapsis. Using the solution of the Kepler problem and 
some constraints in orbital elements, we develop an ex¬ 
tension of the original Hohmann model considering el¬ 
liptical and non-coplanar orbits. We investigate orbital 
changes between elliptical orbits using non-tangential 
impulses which are applied at periapsis and apoapsis 
of the orbit in order to obtain minimum cost of fuel. 
In this paper we find sev eral minimum solutions a nd 
we determine which case ( Kamel and SolimanI 199^ is 
optimal from an energetic point of view. The present 
paper is organized as follows. In Sect. [Hwe discuss the 
problem and the physical constraints used to solve it. 
In Sect. [3] we use a standard optimization technique 
for each configuration and get a polynomial function, 
whose solution (once we found its inverse) gives the an¬ 
gle between the initial orbit and transfer orbit. In Sect. 
[3]we consider some numerical values in order to use our 
solution in a particular example whereas in Sect. Owe 
discuss briefly the solution and consequences. Finally, 
in Sect. Owe present some relevant conclusions of this 
paper. 


2 Problem to solve: A general approach 

In order to visualize all possible trajectories of an or¬ 
biter we consider the initial, transfer, and final orbits, 
with the same nodal line (see Fig. [T|). The initial orbit 
is where first impulse is applied and has a set of orbital 
parameters given by 0{a, e,i,fl,u!,T)i. The second or¬ 
bit, the transfer orbit, is described by a set of param¬ 
eters 0{a,e,i,il,uj,T)t', this orbit is where the second 
impulse will be applied. The arrival orbit has orbital 
parameters 0(a, e, i, H, w, r)y. We want to emphasize 
that each orbit lies in a plane and between any two 
planes it is possible to define an inclination angle. This 
is an important point since we want to find the mini¬ 
mum value given by the minimization of cost function 
(usually the cost function is defined as the sum of im¬ 
pulses per unit mass), taking one of the orbital elements 
as parameter. This will be commented on in Sect. |3l 

2.1 Model Constraints 


• The radius of the initial impulse is lower than the 
radius of the final impulse. 

• The initial, Oi, and final. Of, orbits form an incli¬ 
nation angle iif G (0°,180°), between their orbital 
planes Voi and Vof, respectively. 

• When we apply the first impulse Aui = ||Ax/i||, 
where Aui is the vector of the first maneuver, it pro¬ 
duces an inclination angle in 0°, between the Voi 
and Vot plane of the transfer orbit, O*. 

• When we apply the second impulse Au 2 = ||A'!; 2 ||, 
where Ax /2 is the vector of the second maneuver, it 
produces an inclination angle x*/, between the Vot-, 
and Vof■ 

• The apsides of the three orbits are collineal with the 
nodal line of the three orbits. 

• The primary focus, Fp, is common in the three or¬ 
bits and coincides with the origin of inertial reference 
frame. 

In Fig. [T]we show the three angles iif, in and Uf- Now, 
the angle iif is fixed and defined (because we know in 
advance the angle between initial and arrival orbits), so 
we will choose one of the angles in and itf as the pa¬ 
rameter to minimize. In order to do the minimization, 
we need to define a function that allows us to calculate 
the angles previously mentioned. This is the necessary 
energy for the orbital maneuvers. We call this func¬ 
tion F the cost function] it may depend on all orbital 
parameters. 


2.2 Mathematical aspects 


Hohmann ( I960ll found the transfer of minimum cost 


between two circular orbits using an elliptical transfer 
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Fig. 1 Voi is the plane (blue color) where the initial orbit 
is located, Vot is the plane (green color) that contains the 
transfer orbit. Vof is the plane (red color) that contains the 
final orbit. The arrow directions indicate the motion of the 
satellite in each orbit. Zq. is the so-called line of apsides for 
the initial orbit. Note that it is not required for the initial 
orbit to be interior to the projection final orbit we just take 
r/ 3> )"i as a particular case. 


induces a simple solution, 



Using and the standard impulse definition 

we get the two impulses applied at perigee and apogee, 
and we define the cost function T as the sum of this 
two maneuvers. In the following section we will show 
this procedure in detail. 

3 Method and results 

We define the cost function T as 

^■=11 Aui II + II Au2 II, (4) 

where these impulses (per unit mass) are non-tangen- 
tial, applied at apside extreme lines, that is perigee and 
apogee, respectively (see Fig. [a. The impulsive ma¬ 
neuver vectors are given by 


orbit. In this work we choose the inclination in as pa¬ 
rameter. However, it is possible to choose any orbital 
parameter in order to find a minimu m of cost fun ction 
and get the best possible trajectory (jAbadl 120121) . We 
need to write down the cost function in terms of the 
two impulses. In order to compute the impulses, we 
need to get the norm of each one and relate them with 
the orbit al parameters. This is given by the Vis Viva 
equation (IMontenbruck and Gillll2005ll . 


e 


ii = 1 2 _ ^ 

2a 2 r 


( 1 ) 


where e is the orbital energy, fi = GMq is a constant 
(with Mq Earth mass), r is the relative distance be¬ 
tween the two bodies, and v = ||u||, where v is the 
required velocity vector. Since the two impulses will be 
applied, one at the perigee and the other at the apogee, 
we need to obtain the velocities in perigee and apogee 
for each orbit. Using the formulation for a general conic 
section we get, in polar coordinates, 


ra = (l-|-e)a, (2a) 

Tp = (1 — e)o. (2b) 

Considering Eqs. (1^ . (l2bl) . and © it is possible to 
obtain the norm of velocity in perigee and apogee which 


AFi =Ui-Vi, (5) 

Av 2 =Vf - Uf, ( 6 ) 

where the vectors u) and Ui refer to the initial veloc¬ 
ities of the initial and transfer orbits, respectively. In 
the same way, Vf and Uf are the final velocities of the fi¬ 
nal and transfer orbits, respectively. The first impulse, 
Aui = ||Aui||, is applied in the initial orbit and the 
second impulse is applied in the transfer orbit, neces¬ 
sary to switch to the final orbit. When the first impulse 
is applied, the initial and transfer orbits coincide; that 
happens just in this point, as is seen in Fig. [51 for dif¬ 
ferent configurations. For the first configuration (upper 
left panel of Fig. the initial impulse Avi is located at 
the perigee of the transfer elliptic orbit and this co¬ 
incides with the initial orbit, whereas the final impulse 
Av 2 is located at the apogee of the transfer elliptic orbit 
and this coincides with the final orbit. Note that 
the other three configurations are easily obtained using 
the corresponding impulse maneuvers according to Fig. 
121 Thus, a link between the parameters of the orbits 
is established. A similar situation occurs between the 
transfer and arrival orbits when the second impulse is 
applied. We found an expression for the impulse ma¬ 
neuvers in terms of the velocity vectors for each orbit 
and angles between orbital planes; their norms are 


II AUi II = 

Juf -I- vf - 2uiVi cos(Ut), 

(7) 

II AFz II = 1 

Jvf + uj - 2ufVf cos(U/ - lit), 

(8) 
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Fig. 2 Configurations considering two non-tangential impulses applied in perigee and apogee orbit (only four cases). We 
show constraints relative to this for each case. The figure in the upper left panel corresponds to our first configuration, 
the upper right figure corresponds to the second. The lower left and right figures a re the third a nd fourth configurations, 
respectively. These figures are a plane projection of the three-dimensional problem l|Lacru 312013). 


where Ui = ||Mi||, = Huill, u/ = Hu/H, and Vf = ||u/||. 

Notice that we used the fact iif = in + itf- To find 
the minimum value for the iF function it is necessary 
to calculate the derivative with respect to one of the 
orbital parameters or another free parameter. As we 
already said, we can choose any parameter we want. 
The variation respect to in has been little studied in 
the literature, so we decided to work this case. In order 
to find the minimum value with respect to in for the 
function it is required that 


with ICi and ICf being the kinetic energies before and 
after applying the impulses, per mass unit, respectively. 
To facilitate the algebra we make the following defini¬ 
tions: 


a = cos{in), 

(11a) 

- = sin^(jit). 

(11b) 

in) = Aia + A 2 \/l - cr^. 

(11c) 

in) = Al + + A4a \/1 - 

(lid) 


dF 

din 

F^F 


= 0 , 

> 0 . 


Putting ([T]) and (|i 
a convenient way 

UiVi sin(iit) 


(9a) 

(9b) 

into (I9al) and rearranging terms in 


Fi ‘2uiVi cos(^tn^ 


VfUf sin(zi/ - in) V Ff - 2ufVf cos(zi/ - in) 


, ( 10 ) 


where the Ae set, with .^ = 1, 2,3,4, is given by 

Ai = {cos(fi/), sin(ij/), - cos(2ij/), - sin(2ij/)}. (12) 

Substituting Eqs. (jlla|l . (jllbll . (Illcll and (|lldp . into 
m and aft er some calcul ations we obtain a polynomial 
of the form Lacrud ( 20101) 


j=3 

E 

j=0 


i=6 




j=4 


(13) 
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where the constant bj are 

bo = ujvf{vj + u})- v}u}{u} + vl) cos^(ij/), (14a) 
bi = 2v^u^j:UiVi cos^ {iif) — 2u^v^VfUfCOs{iif), (14b) 

62 = v'fu'j{uf + v1) cos(2iif) — u1v‘l{v'^f + Wj), (14c) 

63 = 2u'lv‘fvfUf cos{iif) — 2vj^uj:UiVi cos{2iif), (14d) 

b 4 = 2ujvfvfUf sm{iif ), (14e) 

b 5 =-vjujiuf+vf)sin{2iif), (14f) 

be = 2v‘fu‘juiVi sin(2zi/) - 2u^vfvfUf (14g) 

Squaring the Eq. (US and regrouping terms, we obtain 
the following sixth-degree polynomial function 


1=6 


1=0 

(15) 

where the Cj are 

Co = - 64, 

(16a) 

Cl = 2(6o&i — bibA ), 

(16b) 

C2 = 2(6o&2 - b ^ bA ) -1-64-1-64-65, 

(16c) 

C3 = 2(6063 — 6560 + 6162 + 6465), 

(16d) 

C4 = 2(6163 -h 6460) + 62 -h 6g — 6g, 

(16e) 

C5 = 2(6263 -l- 6565), 

(16f) 

Co = bl + bl , 

(16g) 


which are represented in an implicit form in terms of 
the velocities and angles, respectively. Now we need to 
obtain the roots of Eq. (HU, since they will provide the 
angle that minimizes the cost function. This procedure 
is similar in each case, but some important differences 
appear due to t he transfer paramet ers eccentricity and 


semimajor axis (IKamel et al.N20111) . In order to keep 


this discussion in general terms, we need to define both 
parameters appropriately: 


(n) 


An) 


i ((1 -h kf\f)af -k (1 -f , 


(1 -I- k^p^ef)af - (1 + k'f'’ei)ai 
(1 -I- k^p^ef)af + (1 + kl"'’ei)ai 


(n) 


(17) 

(18) 


where and are parameters that change depend¬ 
ing of each case (see Table [T] for specific values). In the 
same way, the initial and final velocities are different 
depending on the configurations. All the possible cases 
are shown in Tabled according to Eig. [2] 


4 Numerical example 

In order to illustrate the solution of this problem (i.e. 
find the minimum value of the cost function with re¬ 
spect to the lit variable), we get the cost function for 


Table 1 Selective parameter and for each case. 



Case 1 

Case 2 

Case 3 

Case 4 

IN' 

-1 

-1 

-hi 

-hi 

nW 

-hi 

-1 

-1 

-hi 


Table 2 Initial and final velocity vectors for each case, 
where Upg., and are the velocity vectors in the 

perigee in initial, transfer and final orbits, respectively. 
Vao., Uaoj and Vaoj. are the velocity vectors in the perigee 
in initial, transfer and final orbits, respectively. 


Case 1 Case 2 Case 3 Case 4 


Vi 

^POi 

^POi 




Vaoj 

Vpof 

VpOf 

Vaoj 

Ui 

Vpot 

Vpot 

Vaot 

Vpot 


Vaof 

Vaof 

VpOf 

Vaof 


each case using some known reference values. Eirst of 
all, we need to get two impulse maneuvers to produce 
Hohmann transfer with orbital plane change. For this 
purpose in general six constants are required: /i 0 = 
GMq = 3.98 X 10^^ km^ s“^, the initial-final angle 
iif = 7r/2 rad, and the set {e^, a^, e/, a/} where we 
have taken = 0.052, a,; = 6948 km corresponding to 
the Sputnik I satellite ( NASA|E957 1. and e/ = 0.190, 
a/ = 8682.5 km c orresponding to the Vanguard I satel¬ 
lite (jNASAl 119581 1. Using previous parameters, our so¬ 
lutions give different roots where one of this is the min¬ 
imum global in the cost function between all possible 
roots (six in the most general case). The roots are in 
Table [3] displayed case by case. For instance, in the 


Table 3 Inclination angle obtained by two different tech¬ 
niques: using our sixth-degree polynomial function {iit)p 
and by the numerical solution indicates abso¬ 

lute error in radians. The angles obtained using our proce¬ 
dure and the numerical solution are given in radians. 



Case 1 

Case 2 

Case 3 

Case 4 


0.0577968 

0.0269421 

0.0260117 

0.0505942 

i^it) e 

0.0577970 

0.0269422 

0.0260117 

0.0505981 


0.0000002 

0.0000001 

0 0(7) 

0.0000039 


Case 1 {iit)p = 0.0577968 rad and {iit)e = 0.0577970 
rad, equivalently 3.31152417 degree and 3.3115127 de¬ 
gree, respectively. Therefore, the absolute error is 
A(fit) = \{iit)p - {iit)e\ = 10“'" degree. In Table El 
the term 0(7} means the exact value and the numer¬ 
ical values are in agreement to the seventh order of a 
polynomial expansion of the correct result. 
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5 Discussion 

Since all cases share the same initial and final orbital 
elements, the efficiency between them can be compared. 
In this sense, we note that the cases 1 and 4 (such as 
the cases 2 and 3) can be directly compared, since the 
cost function looks similar. It can be understood due 
to cases 1-4 and 2-3 presenting some symmetry degree 
(reflection). In spite of it, the cost functions of case 1 
and 4 (and case 2 and 3) are different because of the 
particular start and arrival points in the orbit (Fig. [^. 



iit (rad) 

Fig. 3 Cost functions (by case) J- put together as a func¬ 
tion of the initial-transfer angle iu- Black points are the 
minimum values of the cost function for each case. The 
fourth case is the most economical case relative to the other 
three cases. 

We find that two cases are the ones most efficient, the 
first case has the initial impulse at pericenter, with 
the corresponding final impulse at apocenter whereas 
the fourth case has the initial impulse at apocenter, 
with the corresponding final impulse at apocenter. On 
the other hand, cases 2 and 3 are the ones worst with 
respect to the energy used to change between orbits. 
By revisiting the cost function plot we note that the 
more economical configurations arrive at the apocenter 
whereas the more expensive cases arrive at the pericen¬ 
ter (Fig. [3]). Furthermore, as we need to apply an extra 
impulse to stop the vehicle from the transfer orbit to 
final orbit (case 2 and 3), a difference between the cost 
functions is established. 


for the first configuration and calculate the solution 
for the other three possi ble cases. We compared the 
only four possible cases ( Kamel and Soliman I999ll of 
the orbital transfer considering two impulses applied in 
perigee and apogee, and we determine the best model 
for a given set of orbital elements. We compare our 
analytical solution using orbital data from well-known 
satellites. Comparing the exact solution and our solu¬ 
tion, we observe that the roots where the cost function 
is minimal are approximately equal so our result is con¬ 
sidered valid. In the same way, in accordance with our 
solution, the fourth case is the optimal possible case. 
Finally, we show that it is always cheaper to have the 
final impulse at apocenter. 
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6 Conclusions 


In this paper we obtain the analytical minimum cost 
function of an orbital transference between two non- 
coplanar elliptical orbits considering as free parameter 
the inclination between the initial and transfer plane . 
We recovered the solution obtained by Lacruz ( 20ir)ll 
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